Domain walls with non-Abelian orientational moduli 
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Abstract. Domain walls with non-Abelian orientational moduli are constructed in U(N) 
gauge theories coupled to Higgs scalar fields with degenerate masses. The associated global 
symmetry is broken by the domain walls, resulting in the Nambu-Goldstone (and quasi-Nambu- 
Goldstone) bosons, which form the non-Abelian orientational moduli. As walls separate, the 
wave functions of the non-Abelian orientational moduli spread between domain walls. By taking 
the limit of Higgs mass differences to vanish, we clarify the convertion of wall position moduli 
into the non-Abelian orientational moduli. The moduli space metric and its Kahler potential of 
the effective field theory on the domain walls are constructed. We consider two models: a U(l) 
gauge theory with several charged Higgs fields, and a U (N) gauge theory with 2N Higgs fields 
in the fundamental representation. More details are found in our paper[T]. 



1. BPS solitons 

Solitons are useful to build unified models with extra dimensions, and to provide all or part 
of nonperturbative effects. If a global symmetry of the theory is spontaneously broken by the 
presence of solitons, Nambu-Goldstone (NG) bosons come out and form (a part of) the moduli 
space of the soliton. In the case of non-Abelian global symmetry, the resulting massless modes 
can have non-Abelian orientational moduli. Quite often solitons have parameters, which are 
called moduli. If we promote these parameters to fields on the world volume of the soliton, they 
become massless fields in the low-energy effective field theory on the soliton. 

Simplest of solitons is the domain wall which depends only on one spatial dimension, namely 
co-dimension one soliton. In order to have a domain wall, we need to have discrete vacua. As the 
simplest theory with two discrete degenerate vacua, let us consider a 1 + 1 dimensional theory 
of real scalar field (p with a double well potential 



If there is a field configuration connecting the two discrete degenerate vacua, ip + = v and 
= —v, we obtain a domain wall separating two vacua, as a kink, whose energy density is 
localized, resulting in the domain wall. The nontrivial boundary condition at infinity assures 




A > 0. 




the topological stability of the configuration: topological charge is characterized by 7To(A4). 
To obtain such a solution, let us assume a static configuration depending only on one spatial 
direction y, and form the following complete square of the energy density 
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This is called the Bogomol'nyi completion giving the lower bound of energy which is called the 
Bogomol'nyi bound 
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The Bogomol'nyi bound is saturated if and only if the following first order differential equation 
is satisfied: 

d v (p + V\(ip 2 - v 2 ) = 0, (4) 

which is called the Bogomol'nyi-Prasad-Sommerfield (BPS) equation for the domain wall. The 
BPS equation is easily solved to give the BPS domain wall solution 

ip = vtanh(V\v(y - y )). (5) 

The BPS domain wall solution has a single modulus yo, whose physical meaning is the position of 
the wall. This modulus can also be understood as a NG mode corresponding to the spontaneously 
broken translational symmetry. We can promote the real scalar to complex scalar, and add a 
fermion with appropriate interactions to make the system supersymmetric. Namely we can 
embed the theory ([T]) into a supersymmetric theory with four supercharges. The BPS solution 
preserves the half of the supersymmetry in this embeded theory. This is a typical example of the 
BPS solitons which can be understood as the BPS state preserving a part of the supersymmetry 
in a supersymmetric theory. 

As another example of solitons, we can consider vortex in a theory with the Abelian gauge 
field Wfj, coupled to a charged complex scalar field <f> 

c = ~^2 F ^ Fflu + ~ j - y2 Y > ( 6 ) 

where V^cfi = (<9 M + iW^)(p denotes the covariant derivative, and F^ u = dpW v — d u W^ the 
field strength. The mapping from the infinity in the x,y plane (S 1 ) to the vacuum manifold 
\<p\ = v, S 1 gives a topological charge 7Ti(S' 1 ) assuring the stability of the vortex configuration. 

k = -^J d2xF ^ ( 7 ) 

The position of a single vortex is again a modulus, since the vortex can be formed anywhere 
in the x, y plane. If there are two or more vortices, the gauge field induces repulsive force 
between the vortices, and the scalar field induces attractive interaction^]. When e 2 = A, the two 
interactions cancel each other and there is no static force between vortices. Therefore vortices 
can be placed at anywhere relative to each other, and the relative positions become additional 
moduli. This critical case corresponds precisely to the case where the theory can be embedded 
into supersymmetric theory by adding fermions appropriately. The vortex solutions can then 



This applies to vortices with vorticity of the same sign. If the signs of vorticities are opposite, both gauge field 
interactions and scalar interactions become attractive. 



be understood as the BPS solitons preserving the half of the supersymmetry. If there are 
internal global symmetry which are broken by the presence of solitons, the NG modes emerge. 
In particular, if the global symmetry is non-Abelian, we can have non-Abelian orientational 
moduli for the soliton. These non-Abelian orientational moduli exhibits interesting properties 
similar to the D-branes in string theory. 



2. BPS equations for U(N) gauge theories and the moduli matrix 

We consider U(Nq) gauge theory in space-time dimension d = 4 + 1 with a real scalar field 
E in the adjoint representation and N-p (> Nq) flavors of massive Higgs scalar fields in the 
fundamental representation, denoted as an Nq x Np matrix H. Choosing the minimal kinetic 
term, we obtain 
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where the covariant derivatives and field strengths are defined as T)^E = d^E + i[W^,S], 



= (d^ + iWfi)H, F^ v = —i[T>^, T> v \. Our convention for the space-time metric is 
Vnv = diag(+, —,-•-,—). The scalar potential V is given in terms of a diagonal mass matrix M 
and a real parameter c as 
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The 1/2 BPS equations for domain walls interpolating the discrete vacua can be obtained by 
usual Bogomol'nyi completion of the energy 
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The first order differential equations for the configurations saturating this energy bound are of 
the form 17 
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Here we consider static configurations depending only on the y-direction. 

Let us solve these 1/2 BPS equations. Firstly the first equation can be solved by [7] 

H = S-\y)H e M y, E + iW y = S- l (y)d y S{y). 
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Here Hq, called the moduli matrix, is an Nc x Np constant complex matrix of rank Nq, and 
contains all the moduli parameters of solutions. The matrix valued quantity S(y ) £ GL(Nq, C) 
is determined by the second equation in (|13|) which can be converted to the following equation 
for U = SS^: 
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This equation is called the master equation for domain walls. From the vacuum conditions at 
spatial infinities y — > ±oo, we can see that the solution 17 of the master equation should satisfy 
the boundary condition $7 — > Qq as y — > ±oo. It determines S for a given moduli matrix Hq 
up to the gauge transformations S~ l — ► US~ 1 , U G U(Nq) and then the physical fields can 



be obtained through (|14p . Note that the master equation is symmetric under the following 
^-transformations 

Hq^VHq and S(y) -> VS(y) with V G GL{Nq, C), (16) 

and if the moduli matrices are related by the ^/-transformations Hq = VHq, they give physically 
equivalent configurations. We call this equivalence relation as the ^-equivalence relation and 
denote it as Hq ~ VHq. The solution of the master equation exists and unique for any given Hq 
at least for the U(l) gauge theory [5]. For U(N) gauge theory the number of moduli parameters 
agrees with the the result of index theorem [6]. 

3. Non-Abelian orientational moduli of walls in a U{1) gauge theory 

So far, domain walls with eight supercharges have been mostly considered in gauge theories with 
U{1) gauge field [SJ-JSJ, U(l) x U(l) gauge field [II], or U(N) gauge fields [7]-[IU] coupled to 
Higgs scalar fields with non- degenerate masses except for [12|. I13j. In the case of non-degenerate 
Higgs masses, the flavor symmetry is Abelian: U(1) Nf ~ 1 and the symmetry of the vacua is also 
Abelian. As a result each domain wall carries a U(l) orientational modulus [21 [7]; The moduli 
space of a single domain wall is 

W fe=1 ~ R x S 1 . (17) 

From this viewpoint, these domain walls should be called Abelian domain walls even when the 
gauge symmetry of the Lagrangian is non- Abelian [7]-|llj. 

Let us see the non-Abelian orientational moduli in a simple example of the Abelian gauge 
theory coupled with the Np = 4 Higgs fields. 

3.1. Vacua of the U(l) gauge theory 

The massless vacuum manifold is T*CP 3 where the base manifold is parametrized by 

CP 3 = [HH* = c} /[/(l), H = yfc(h l , h 2 , h 3 , h 4 ) , (18) 

where the quotient is the overall U(l). The vacuum manifold is expressed as (the inside and 
the surface of) a triangular pyramid in the 3 dimensional space (|/ii| 2 , |/i2| 2 , |/i3| 3 ), as shown in 
Fig- H (c). When the mass matrix containing a small parameter e (0 < e G R) 

/ me me \ 
M = diag Im, — , — — , -mj (19) 

is turned on, the vacuum manifold is lifted except for four points and the flavor symmetry breaks 
from SU(4) to U (l) 3 . These discrete vacua are the four vertices of the pyramid shown in Fig. [H 
(a). We label those vacua as (A) (A = 1,2,3,4). The vacuum expectation value (VEV) of the 
vacuum (A) is /is = 5ab- Taking a limit of e — » 0, the second and the third Higgs fields become 
degenerate so that the flavor symmetry enhances from £7(1) 3 to U(l) 2 x SU(2) G S'C r (4). There 
are two isolated vacua and one degenerate vacuum CP 1 ~ SU{2)/U{1) represented by a line 
connecting (2) and (3) as shown by a thick line in FigQJb). We denote this degenerate vacuum 
as (2-3). 

3.2. Domain walls in the U(l) gauge theory 

There exist domain wall solutions interpolating vacua in the model with fully or partially non- 
degenerate Higgs masses. In the case of Ac = 1, the moduli matrix and the ^-equivalence (|16|) 
take the form of 



Hq = (0i, 4>2, 4>3, 4><i) ~ A (0i, 4> 2 , 4>3, 4><i) , AG C*. 



(20) 




Figure 1. Vacua for various cases of mass configurations plotted in the three-dimensional space of 
Higgs fields hf,i = 1,2,3 with J2i=i = 1- ( a ) non-degenerate massive vacua (b) massive degenerate 
and non-degenerate vacua (c) massless vacuum. 



In terms of the moduli matrix the vacua (^4) are described by 4>b = $BA for B = 1,2, 3, 4. Since 
we want to consider the domain wall interpolating the vacua (1) and (4) (passing by (2) , (3) 
on the way), the parameter 4>i and ^4 should not be zero while cj)2 , ^3 can become zero. So the 
moduli space corresponding to the multiple domain walls which connect (1) and (4) is 

M ~ (c 2 x (C*) 2 ) //C* ~ C* x C 2 , (21) 

where double slash denotes identification by the ^-transformation. Here the part C* ~ Rx(7(l) 
represents the translational modulus and the associated phase modulus. 

When we take the gauge coupling g to infinity, the model reduces to a nonlinear sigma model 
whose target space is the Higgs branch of the vacua T*CP 3 in the original theory. To make the 
discussion simple, we take this limit for a while. One benefit to consider the nonlinear sigma 
model is that the BPS equations can be analytically solved. In fact the solutions are expressed 
as [8] 

H = -^H e My with ft = H e 2My Hl (22) 
v^o 

A domain wall solution corresponds to a trajectory connecting the vertex (1) and (4). Flows 
from (1) to (4) inside the pyramid are shown in Fig. [2j 

Physical meaning of the moduli parameters becomes much clearer by using the 1^-equivalence 
relation (|20p to fix the form of the moduli matrix as 

H = (1, e Vl+V2 , e ^+^+^3) _ (23) 

Furthermore, one may be visually able to see the "kink" configuration in the profile of the field 
S = (l/2)dy log Qq. In the vacuum region (A) the function takes the value S = ijia- 

Several solutions are shown in Fig. [3l The domain wall positions can be roughly read from the 
moduli matrix in Eq. (|23p as 

y~\- 1 * . * y~ 1 * inA\ 

— = <p l + ( Pl , =(f2 + (f2j =(fi3 + (p3i (24) 

L + Lq L- 

where y+ is the position of the right wall and yo, y_ are the positions for the middle and the left 
walls, respectively. Here £±,0 stands for the width of each wall 
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(a) non-degenerate mass (b) degenerate mass 

Figure 2. Domain wall trajectories in the target space CP 3 for non-degenerate mass (a) and for 
degenerate mass (b). 




(a) e = l/3 (b)e=l/20 (c) £ = 1/100 



Figure 3. Configuration of £ (first row) and density of the Kahler metric of ipi, ip2 and ip3 (second 
row). Moduli parameters are (tp±, if2, f3) — (20, 0, —20) and m = 1. 

This rough estimation is, of course, valid only for well separated walls whose positions are 
aligned as y_ -C yo <C y+, see Fig. [3] (a). Each domain wall is accompanied by a complex moduli 
parameter ipi whose real part is related to the wall position and imaginary part is the U(l) 
internal symmetry (the NG mode associated with the broken U(l) flavor symmetry). 

To argue symmetry aspects of the moduli parameters, first let us consider a model which has 
completely non-degenerate masses and domain walls interpolating between those vacua. The 
global symmetry explicitly breaks from SU(A) to U(l) 3 C SU(4). We take, as the unbroken 
global symmetries, ?7i (1), U%(1) and £73(1) with generators diag (1, — 1, — 1, 1), diag (1, 0, 0, — 1), 
and diag (0, 1, —1, 0) respectively. Each vacuum (A) preserves all of these symmetries. However, 
once domain walls connecting those vacua appear, they break all or a part of these symmetries. 
For example, the moduli matrix Hq = (1, 0, 0, (£4) corresponding to a domain wall connecting two 
vacua (1) and (4) breaks J7 2 (l) but still preserves £7i(l) and ^3(1). Here note that overall phase 
can be absorbed by the ^/-transformation (|20p . Therefore the phase of the moduli parameter 
eft 4 corresponds to nothing but the broken global symmetry L^j(l). This implies the NG mode 
localizes around the domain wall as we saw above. For the moduli matrix Hq = (1, (ft%, 0, ^4), 
which corresponds to two domain walls connecting three vacua (1) — > (2) — * (4), the symmetry 
1/3(1) in addition to ^(l) breaks while a combination of Ui(l) and ^3(1) is still preserved. 
Moreover, when we turn on the third element in the moduli matrix as Hq = (1, 02, ^3, ^4), the 
third vacuum region appears and then the configuration has three domain walls connecting four 



vacua (1) — > (2) — > (3) — > (4). In this case all of £/(l) 3 are broken by the domain walls, so that 
corresponding three NG modes appear. These three NG modes are described by imaginary parts 
of log0, which are combined with the three positions (|24p . to form three complex coordinates 
of the moduli space C 2 x C*. 

Next we consider a limit where the second and the third masses are degenerate (e — > in 
the mass matrix (|19p ). In the limit the global symmetry J7i(l) x 1/2(1) x Us(l) is enhanced 
to Ui(l) x 1/2(1) x SU(2). At the same time, the degenerate vacuum (2-3) appear instead 
of the two isolated vacua (2) and (3) as shown in Fig. [T] (b). At the degenerate vacuum, 
U±(l), 1/2(1) are preserved but SU(2) is broken to ^3(1). Therefore the degenerate vacuum 
(2—3) is SU (2)/ 1/3(1) = CP 1 . Non-vanishing 04 7^ causes the wall interpolating two vacua 
(1) — ► (4) and breaks only 1/2(1) again. Once the degenerate vacuum appears in the configuration 
such as two domain walls connecting vacua like (1) — ► (2—3) — > (4), the breaking pattern of 
the global symmetry becomes different from that in the case of fully non-degenerate masses. 
The moduli matrix Hq = (1, 02)03)04) describes such domain walls. Note that the second and 
the third elements break SU(2) completely. The global symmetry U\(l) x 1/2(1) x SU(2) are 
broken to U(l) which is a mixture of Ui(l) and H £ SU(2). Emergence of the second wall and 
further [/(l)-symmetry breaking are related to the facts that |02| 2 + 103| 2 7^ and 04 7^ 0. These 
facts imply that the modes corresponding to the two broken C7(l)'s localize around the walls 
accompanied by the two position moduli and the mode corresponding to SU(2)/H have support 
in a region around the degenerate vacuum (2—3). We can count the number of the moduli 
parameters as follows. Two real parameters {|02| 2 + |03| 2 5 104| 2 } correspond to the positions 
of the two walls whereas remaining four parameters correspond to the broken global symmetry 
Ut(l) x 1/2(1) x SU(2)/U(1). This is again consistent with dim R (C 2 x C*). 

In the Fig. [3] we showed domain wall configurations of the three domain walls connecting 
the four vacua. As the parameter e decreases, the width of the middle domain wall connecting 
the vacua (2) and (3) becomes broad and the tension of the wall becomes small since they are 
proportional to 1/e and e, respectively. When the width of the middle wall becomes larger 
than the separation of two outside walls, Lq > y + — y_, we can no longer see the middle wall. 
The density of the Kahler metric for the moduli parameters tpi, cp2 and in the strong gauge 
coupling limit are shown in the second row of Fig. [3j The Kahler potential in the strong coupling 
limit is given by K = c J dy log [T3]. When three walls are well isolated as Fig. [3] (a), three 
modes corresponding to the moduli parameters tpi, (f2 and if 3 are localized on the respective 
domain walls. As e decreases, the density of the Kahler metric of ip% is no longer localized but 
is stretched between two outside domain walls. In the limit where e — > the physical meaning 
of (f2 as the position and the internal phase associated with the middle domain wall should be 
completely discarded. Instead, (f2 gives the non-Abelian orientational moduli which comes from 
the CP 1 zero modes of the degenerate vacua (2—3). For fixed moduli parameters ^1,^2,^3, 
the domain wall solution as a function of y sweeps out a trajectory in the target space CP 3 . 
These domain wall trajectories are shown for various values of moduli parameters in Fig. [5J the 
non-degenerate mass case (a) and the degenerate mass case (b). For the degenerate mass case, 
the trajectories do not go out from the triangular plane whose vertices are (1) , (4) and one point 
on the edge between (2) and (3). 

4. The Generalized Shifman-Yung (GSY) Model 

4-1- U(N) gauge theory 

Let us now consider non-Abelian gauge theory with degenerate masses of the Higgs fields. 
Previously considered model is the U (2) gauge theory with four Higgs fields in the fundamental 
representation with the mass matrix M = diag(m, m, — m, — m) \12\ I13j . which we call the 
Shifman-Yung Model. The model enjoys a flavor symmetry SU (2)i, x SU(2)n x U(1)a- It has 
been demonstrated that the coincident domain wall configurations break the flavor symmetry 



to SU(2) Y and the NG bosons corresponding to [SU(2) L x SU(2) R x U(l) A ]/SU(2)v ^ U(2) 
appear in the effective action on the walls. 

By generalizing the Shifman-Yung model, we consider the U (N) gauge theory with Np = 2N 
Higgs fields in the fundamental representation whose mass matrix is given by 

N N 

a 3 „ 1 




M = m— ®1tv = -diag(m, • • • , m, — m, ■ ■ ■ , — m). (26) 

This system has a non-Abelian flavor symmetry SU(N)i J x SU(N) R x U(l)\. Since we have 
only two mass parameters m and —m, possible vacua are classified by an integer < k < N: in 
the k-th vacua, there is a configuration in which k flavors of the first half and N — k flavors of 
the latter half take non-vanishing values and then £ and H are 

N-k 

1 * . 

I vacuum 

m /? ( ^ io7\ 

-"I vacuum _ \ N _ k l N _ k )' 

This vacuum is labeled as (k,N — k). The flavor symmetry SU(N)l is broken down to 
SU(k) c +L x SU(N - k) L x C/(l) c +L, and SU(N) R is broken down to SU(k) R x SU(N - k) c+ R x 
C7(1)c+r, where C + L,C + R denote the locking of the L, R flavor symmetry with the U(N)c 
color symmetry. Consequently the number of the discrete components of the vacua is N + 1 in 
this system. Therefore in this vacua there emerge 4k(N — k) NG modes, which parametrize the 
direct product of two Grassmann manifolds, 

Walls are obtained by interpolating between a vacuum at y = — oo and another vacuum 
at y = oo. The boundary conditions at both infinities define topological sectors. For a given 
topological sector, we may find several walls. The maximal number of walls in this system is N, 
which are obtained for the following maximal topological sector 

/ ^/c(1tv,0tv) at y = +00 , > 



^/c(0jv,l7v) at y = -00 

The unbroken symmetries of the vacua (N, 0) and (0, N) which we consider as the boundary 
condition of domain walls here, are SU(N) c +l x ^(1)c+l x SU(N) r and SU(N) L x SU(N) c+ r x 
C/(1)c+r, respectively. In this case, the moduli matrix Hq can be set into the following form 
without loss of generality: 

^o = v^(ljv,e*)~ v^(e- , liv), (30) 

where is an element of GL(N, C) and <fi describes the moduli space of walls of this system, 
and the two forms are related by the ^/-transformation (|16p . Therefore the moduli space of 
domain walls in the GSY model is 

M ~ GL(N, C)[= U(N) C ] ~ [C* x SL(N, C)] /Z N . (31) 

This moduli space admits the isometry 

- e ia g L e*g R (32) 

with (s'LjA'r) S SU(N)l x SU(N) r and e ta G £7(1)a- This is because the domain wall solutions 
break the symmetry of the two vacua (N,0) and (0,N), G = SU(N) c+ l X 5C/(iV)c+R x 
C7(1)c+l-Rj down to its subgroup. 



4-2. Nambu-Goldstone (NG) modes and quasi-NG modes 

Note the fact that the global symmetry G = SU(N) L x SU(N) R x U(1) A in <[32j) acts on 
the moduli space metric as an isometry whereas the complexified group G c = SL(N, C)l x 
SL(N, C)r x C* acts on it transitively but not as an isometry. Therefore G c action may change 
the point in moduli space to another with a different symmetry structure. Since the symmetry 
of Lagrangian is G but not G c we can use only G when we discuss the symmetry structure at 
each point in moduli space. General 4> can be transformed by G to 

= diag.(t»i,u 2 , ■ • • ,v N ) (33) 

with Vi real parameters. When all v^s coincide (coincident walls), the unbroken symmetry is 
the maximal H mSLX = SU(N)y, so we call it the symmetric point. There are iV 2 massless NG 
bosons and iV 2 quasi-NG bosons. When all v^s are different from each other (separated walls), 
iJ max = SU(N)y is further broken down to H m i n = C/(l)y _1 . Here the numbers of NG bosons 
and quasi-NG bosons are 2N 2 — (N — 1) and N — 1, respectively. These N — 1 quasi-NG bosons 
correspond to the N — 1 parameters Vi without the overall factor. Therefore some of quasi-NG 
bosons at the symmetric point in the moduli space change to the NG bosons parametrizing 
H max / H m { n = SU(N)y/U(l)y _1 reflecting this further symmetry breaking. When some fj's 
coincide, some non-Abelian groups are recovered: H = U (l) v x FJ {/(nj)v- Then the NG modes 
iJ max /i? are supplied from quasi-NG modes. 

The diagonal moduli parameters V{ (quasi-NG bosons) in Eq. (|33|) correspond to the positions 
of N domain walls. When all domain walls are separated, the unbroken symmetry is C/(l)y _1 . 
When positions of n domain walls coincide, U(n)y symmetry is recovered. This phenomenon 
has a resemblance to the case of D-branes. However, there is a crucial difference: the symmetry 
in our case of domain walls is a global symmetry, whereas that of D-branes is a local gauge 
symmetry. However in the case of the d = 2 + 1 wall world-volume, massless scalars can be 
dualized to gauge fields. Shifman and Yung [12] expected that the off-diagonal gauge bosons of 
U(N) (which are originally the off-diagonal NG bosons of U(N) before taking a duality) will 
become massive when domain walls are separated, in order to interpret domain walls as D- 
branes. However, our analysis shows that the off-diagonal NG bosons of U(N) remain massless, 
and instead some of the quasi-NG bosons become NG bosons for further symmetry breaking 
with the total number of massless bosons unchanged as explained. 

We have obtained the low-energy effective Lagrangian explicitly for this model. We found the 
following. The wave functions of the NG boson for translation, for U(l) , and quasi-NG bosons 
are localized, and that other massless modes are extended between two domain walls, if walls are 
separated. Wave functions of all the massless modes become identical in the limit of coincident 
walls. These behaviors are different from D-branes, although there are some similarlities. More 
precise details of our results can be found in Ref.pQ. 

When we introduce complex masses for Higgs fields, domain wall junction or network emerge 
as 1/4 BPS states |15|-|17j. In this case too, non-Abelain NG modes appear in the effective 
action, when some masses are degenerate [18], [19]. 
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Appendix 

We will briefly give the derivation of our master equation (|15p . Let us first note an identity for 
variations Si, i = 1,2 of arbitary regular matrix A 

5 1 (A- 1 6 2 A) = A^S&A - A- 1 5 1 AA- 1 5 2 A = A~ 1 (5 2 S 1 AA' 1 - 5\AA~ l 5 2 AA~ 1 )A 



= A- 1 5 2 {S 1 AA- 1 )A 



(A.l) 



For gauge invariant quantity O = SS^ , we obtain 

^(fi&fi -1 ) = 5 1 (ss^5 2 (s-^s- 1 )) = s^ss^i^s-^s- 1 + S52S- 1 ) 

= S (<Ji(5 t «5 2 5- lt ) + S- 1 5 1 SS ] 5 2 S-^ + 6 2 S~ lt SiS' 1 S + S- 1 6 1 (S5 2 S- 1 )S^ S' 1 
= S (^(S^S -1 *) - 6iS- 1 SS^d 2 S-^ + S^6 2 S-^5 1 S- 1 S + d^SxS^S)) S' 1 (A.2) 

By choosing both 5\ and 5 2 to be dirivative in y, and noting that 

£ = + dytfS- 1 *), Wy = ^{S- l dyS - dyS^S^), (A.3) 

we obtain the master equation (|15p for the domain walls. Choosing 5\ and £ 2 as d z and 5^ 
respectively, we can obtain the master equation for vortex. 
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